Pure nuclear quadrupole resonance determination of the electric field gradient 
asymmetry for broad lines (and application to YBa 2 Cu 3 07). 
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We present an angle dependent nuclear quadrupole resonance (ADNQR) method to determine 
the electric field gradient asymmetry parameter rj in systems where the resonance line is so broad 
that the radio frequency field can excite only a portion of the nuclear spins. In this situation the 
recently developed spectroscopic methods are not applicable. ADNQR is useful for single crystals 
and oriented powders, and, for small n determines rj 4 . Therefore, it can be used to evaluate fluctu- 
ations in rj due to inhomogeneities. We demonstrate the application of ADNQR experimentally to 
oriented superconducting YBa2Cu307 powder. 
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INTRODUCTION 



Pure nuclear quadrupole resonance (NQR) is a very 
useful tool for studying electronic properties of materi- 
als without the need to apply an external magnetic field. 
One famous example is in the study of superconductors 
where magnetic fields penetrate the sample in a non uni- 
form way or do not penetrate at all. The NQR is de- 
termined by two parameters: a frequency scale v q and 
asymmetry parameter r\ each carry important informa- 
tion on the charge distribution in the system under in- 
vestigation. However, for nuclei with I = 3/2 there is 
only one resonance frequency out of which one cannot 
separate v q and rj. In recent years some ingenious ex- 
perimental methods have been invented to extract these 
parameters without applying a field. This include ID 1], 
and amplitude Q and phase modulated 2D nutation 
spectroscopy. However, these methods are based on the 
assumptions that the NQR line is much narrower than 
the effective nutation frequency (which is on the order of 
7 Hi where Hi is the radio frequency (RF) field and 7 
is the nuclear gyromagnetic ratio). Naturally, there are 
many occasions where these requirements are not met; 
superconductors, magnetic materials, and systems with 
disorder are just a few examples. 

In this paper we provide an alternative method, based 
purely on NQR, of determining rj (and therefore v q ) sep- 
arately for a spin 3/2 in systems with broad lines. In 
addition, in cases with tetragonal symmetry (77 = 0) , our 
new method is sensitive to the fluctuations in rj, thus pro- 
viding a measure of the system's homogeneity. However, 
this method is applicable only to single crystals or ori- 
ented powders where all observed nuclei have the same 
principal axis of the EFG. The main idea of the new 
method is to measure the signal intensity as a function 
of the angle between the RF field Hi(t) and the principal 
axes of the EFG, hence it is named angle depended NQR 
(ADNQR). In the experiment the sample is rotated with 
respect to the symmetry axis of the coil. 
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FIG. 1: The expected echo intensity in a single crystal as 
a function of the polar angle 8 for a fixed azimuthal angle 
ip. Two extreme values of rj and <p are examined. The pulse 
length is optimized at 6 = 90 and kept constant throughout 
the sample rotation. 



THEORY 

A full theory of the signal intensity as a function of 
RF field direction was given by Pratt 4]. However, his 
focus was mainly on the line shape moments in the case 
of dipolar coupling, and not on the extraction of rj. We 
therefore re-derive the same calculation from our view- 
point. The starting point is the quadrupole Hamiltonian 
given by 
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H q ='-^[3I 2 z -I 2 + r,(I 2 x -I*)] 



(1) 



As mentioned before, its energy levels are two doublets 
with angular frequency separation 



to 
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We further assume that there are no interactions among 
the nuclear spins and that the line width is due to local 
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variations in the values of v q and 77, namely, inhomo- 
geneous broadening. In general, such broadening could 
stamp from crystal defects, impurities, or grain bound- 
aries. In the cuprate superconductors (which will be ex- 
amined here) the broadening might be due to charge seg- 
regation 0. In these cases of broad lines, a potential use 
of ADNQR is to scan ui point by point, determine 77 at 
each point by ADNQR, and extract v q at each point from 
Eq.0 

The RF field is transmitted in the r direction on reso- 
nance only with a portion of the entire line. The expec- 
tation values of the excited spins at time t after the RF 
application is obtained by solving 
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= [Tt q — tkviTr cos(wi)] 



where u>i — jHi . This is done using the matrix 
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where a = rj/[y/3(l 
and gives 
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where the superscript A on any operator O means A = 
A^OA. By introducing a transformation equivalent to 
moving into the rotating reference frame in NMR 



\n) = Ae~ m >' h \ri) , 
and substituting Eq. [S] in Eq. [3] one finds 
d 



(5) 



i lk = ~ UllA ^ ' rcos(wi) \n') 



where 



l A (t) = e< l l h \ A 



,-mft/n 



Keeping only the secular (time independent) terms of 
I A (t) and calling them I A , an operation equivalent to 
ignoring the terms contra-rotating at a frequency 2uj in 
NMR, leads to the simplified equation 



and its solution 



— Loil A r \n') 
exp{jw\l A • rt) \n') . 



The RF field is applied for a time period t^/ 2l after which 
the spins evolve according to Tl q only. Therefore 



AA^e- m " {t -^^ )/h Ae 



exp(iojil A -rt n 



namely, 



where t is the time measured from the moment the pulse 
started. In a single coil FID experiment the signal at 
time t is given by 



(Lr) 
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Tr exp(-/3Wf ) expH^P 4 ■ rt n/2 ) (6) 



[I A (t) •r]exp(iwiI A -rt n /2), 



where Z is the partition function, [3 = l/fc^T, and T is 
the temperature. However, here again we are interested 
only in terms oscillating at frequency uj. This means that 
we can replace I A (t) by 

I A (t) - F(i) = ^Kt/n^A e -iHft/n_ 

I A (t) contains only the oscillating terms at frequency u> 
inI A (i). 

For a general irradiation direction, Eq. gives in the 
high T approximation 
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(Lr) t = — - sin(Awi^ /2 ) sm(wrj) (7) 



where A is an efficiency factor given by 



with 



(77 + 3,77-3,277). 



(8) 



For an echo experiment, obtained by a7r/2 — t — tt pulse 
sequence, the phase accumulated during the 7r/2 affects 
the echo intensity as the 3rd power of its sinus, as demon- 
strated for completion in appendix B This leads to our 
major finding 
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sin (Awi^/ 2 ). 



(9) 



From this point on we will be interested only in 
the echo amplitude A at time 2r as a function of 
the polar and azimutal angles 9 and ip where r = 
(sin 9 cos ip, sin 9 sin ip, cos 9) . We therefore define 



A(0,<p) = {I# 



Echo 



There are two particularly interesting cases: tp = (xz 
plane) and ip = 90 (yz plane), for which 



/2J 



\(9,<p = Q) = 



(77 + 3) 2 sin 2 9 + 4r] 2 cos 2 9 



4(3 + rf) 
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and 



X(9,cp = 90) = 



'(7?-3) 2 sin 2 6» + 4?7 2 cos 2 9 
4(3 + 77 2 ) " 



The echo intensity as a function of 9, for <p — and 
tp = 90 and for rj = and 1, normalized by the intensity 
when = 90, is depicted in Fig. ^ The n/2 pulse length 
is optimized at 9 = 90 (by \w\t v /2 = ir/2) and then kept 
constant throughout the sample rotation. For 77 <§C 1 the 
echo disappears when the RF is applied in the z direction 
(9 = 0). In addition, there is no difference between the 
9 dependence of the echo in the two planes, ip = and 
ip — 90. On the other hand for rj = 1 there is different 
behavior in the two planes. Along the yz plane ((p = 90) 
the echo intensity is constant. In contrast, in the xz plane 
{ip = 0) it drops to 27.6% of its maximum strength. 

An important limit which could be further analyzed 
is 1] <C 1. Since there is no difference between the 
two planes we denote the efficiency factor at 9 = 90 
by Aj^ = This result was first obtained by Das 

and Hahn [6|. In this case the optimum pulse length is 
given by t w /2 = 7r /u>iy3 at which the signal amplitude is 
A(90, p) — (3uj^/3/A. When keeping the pulse length con- 
stant and transmitting in any other direction the relative 
signal intensity will be 

MM = 2A . 3 ttA 
A{99^) V3 &ml V3 j ' 

In particular, when transmitting along the z direction 
A = a z ~ |ry| /\/3, and 



A(90,p) 
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This result is important for cases where r\ fluctuates due 
to crystal inhomogeneity. To appreciate such fluctuations 
one can determine the averaged 77 by standard methods 
and the fourth moment of the r\ distribution by ADNQR. 
Thus measurements of A(0, p) could provide information 
on the quality of the crystal. 

Next we discuss the situation of oriented powders 
where the z direction is well defined but the xy direc- 
tions are interchangeable. One example is high temper- 
ature superconductors. The samples are usually made 
from grains of single crystal oriented in a magnetic field 
so that their z direction aligns with the field direction. 
However, the x and y directions are random. In this case 
the outcome of our experiment is a planar average of the 
result presented in Fig. ^ namely, 

2tt 



A{9) = : kl A ^>* 



This averaging can only be done numerically. In Fig.[21we 
present the angular dependence of the echo intensity in 
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FIG. 2: The expected echo intensity in an oriented powder as 
a function of the polar angle 9 for various values of y. The 
azimuthal angle ip is averaged. The pulse length is optimized 
at 6 — 90 and kept constant throughout the coil rotation. 
The inset shows the ratio of the echo intensity between 6 = 
and 9 — 90 as a function of y. 



the case of oriented powder. Here again the pulse length 
is chosen so that A(90) is optimized. In the i] = 1 case, 
37% of the signal remains when Hi is in the z direction 
compared to the signal when H\ is in the xy plane. In 
contrast, in the 77 = case the signal is lost completely 
under the same experimental conditions. In the inset of 
Fig. [3 we show A(0)/A(90) as a function of 77. It is clear 
that A(0) is most sensitive to rj > 0.5. 



EXPERIMENT 

We demonstrate an experimental application of this 
method to 63 Cu in oriented powders of the high temper- 
ature superconductor YBa2Cu307 at room temperature. 
This system has an orthorhombic symmetry, yet all the 
reports using NMR are consistent with 77 = 8J. The 
powder samples arc mixed in a teflon container with sty- 
cast 1260 and placed in a magnetic field of 9T at room 
temperature for 24 hr. During the first 1/2 hr it is shaken 
with a motor. After drying, a solid cylinder 1.0 cm in di- 
ameter and a 2.0 cm long sample is removed from the 
teflon container. In order to obtain the NQR line, the 
sample is placed in a coil which is tightly wound around 
it. The NQR line is measured with a Tecmag Apollo 
spectrometer to which a home-built automated frequency 
sweep feature is added; the circuit remains tuned and 
matched throughout the frequency sweep. We find that 
the NQR line is centered at 31.6MHz and has a 1MHz 
width. These properties and the line shape are in agree- 
ment with those previously published 0- In addition, T2 
is on the order of 50 fi sec and depends on the frequency. 
For the rotation measurements the sample was placed 
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ECHO INTENSITY FOR NON PERFECT PULSES 
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FIG. 3: Experimental demonstration of the variation in 
echo intensity as a function of 8 in oriented powders of the 
YBa2Cu307 superconductor. The solid line is the theoretical 
expectation for r\ = 0. 



In order to appreciate the influence of a non perfect 
pulses on the echo intensity, we examine the situation for 
a single spin 1/2. In a two pulse experiment, in which the 
second pulse is twice as long as the first one, the nuclear 
spin expectation value at the time of the echo is given by 



(^> 2r = Tr{ 



exp(-PH) 



cxp(—il x 9) exp(—iLuI z T) exp>(—2il x 9) exp>(—iujI z T)I x 
exp(iwJ z T) exp(2il x 6) exp(iu)I z T) exp(il x 9)} 

where 9 — I n the high temperature approxima- 

tion this simplifies to 



(7 x ) 2r = ^ s\n9[s\n 2 9-cos 2 9 {2cos{uj t) + cos(w 2r)}] 



across a 5 cm long coil to achieve better field homogene- 
ity, and was rotated by hand. In these measurments only 
5% of the line at its center is excited. 

We present the rotation data in Fig. [3] We varied 9 
over 360 degrees to demonstrate the fact that A(9)/A(90) 
is always positive and periodic. The signal intensity 
varies by a factor 20 between the 9 = and 9 = 90 
directions. However, our numerical calculations of Fig. [2 
(solid lines) do not fit the experimental results exactly. 
There could be three possible reasons for this. One is 
that the RF field inside the coil is not uniform enough, 
the second is that YBa 2 Cu307 is not perfectly oriented, 
and the third is that, while the average rj in YBa2Cu307 
is zero, fluctuation are important due to charge inho- 
mogeneity. Further study is required to address these 
possibilities. 



CONCLUDING REMARKS 

We have demonstrated that ADNQR could be used 
to determine rj and therefore v q at particular frequen- 
cies along a broad NQR line. In a future publication we 
will use ADNQR to determine the charge fluctuations 
evolution as a function of doping in high temperature 
superconductors. 
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The oscillating term which depends on r is a residue of 
the second pulse, and the one which depends on 2r is a 
result of the first pulse. They contribute only if these 
pulses are not perfect, namely 9 ^ 7r/2. In reality, these 
oscillating terms will relax before the echo is formed even 
if the pulses are not perfect. So if r is long enough 
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This is the origin of the cubic pulse length dependence 
of the signal as a function of the RF field. 
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